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On weak convergence of stochastic heat equation with 

colored noise* 

Pavel Bezdek^ 


Abstract 

In this work we are going to show weak convergence of a probability measure 
corresponding to the solution of the following nonlinear stochastic heat equation 
^ut{x) = ^/!\ut{x)+a{ut{x))'qa with colored noise rja to the measure corresponding 
to the solution of the same equation but with white noise ry as a t 1 on the space of 
continuous functions with compact support. The noise is assumed to be colored 
in space and its covariance is given by E [ria{t, x)r]a{s, y)] = 6{t — s)fa{x — y) where 
fa is the Riesz kernel fa{x) oc 1/ |x|“. We will also state a result about continuity 
of measure in a, for a G (0,1). We will work with the classical notion of weak 
convergence of measures. 

Keywords : The Stochastic heat equation, colored noise, Riesz kernel 


1 Introduction 

Throughout this work we will consider the following one-dimensional heat equation 

Q Kj 

-^Ua,t{x) = -k (T(Ua,t{x))r]a , X G R , f > 0, (1) 

with colored noise having a particular covarinace structure 

E x)ry„(s, y)] = 6{t - s)/«(x - yy), (2) 
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where 


fa{x) = ci_a.ga{x) 


1 


for a G (0,1), 


( 3 ) 


= gi-a[x) , go,[x 


a 


X 


and the constant Cq, is 



( 4 ) 


The initial condition Ua,o(x) is taken to be bounded and measurable. We will also assume 


a to be Lipschitz continuous with Lipschitz constant K. Stochastic PDEs such as ([1]) 
have been studied in [71 [I2l [21 lUl E] and others. 

The function /„ can be thought of as an ‘approximation’ to the delta function in 
the following special sense, we know that one-dimensional Fourier transform of gi-a , 
denoted by gi-a, is equal to fa- We also know that the Fourier transform of a constant 
is 6 distribution. Observe that gi-a converges pointwise to 1 as a 11- We will study the 
solution of dl]) as a function of a. This arises noticeably in [H Sec. 7] where the authors 
have shown that Lf{P) norm of Ua,t{x) converges to Lf{P) norm of the solution to ([5]) 
as a 7 1 for every t > 0,x eM and a{x) = x. 

The main question that has motivated this work, is whether the solution of ffTl) 
converges [in the appropriate sense] to the solution of the same equation, but with white 
noise g instead of colored noise as a 7 1- By that we mean, the solution to 



( 5 ) 

( 6 ) 


uq{x) = Uaft bounded and measurable. 


where g denotes white noise. We will state the main theorem in terms of measures 
corresponding to solutions. Let C = C([0,T] x [—A^, iV]) for any hxed compact set 
[0,T] X [—iV, iV] C Mq X M. Denote Pa, the measure corresponding to Ua, 



P {ua G A} for a G (0,1), 

P {m G A} for a = 1, 


for any Borel set A of space C. Here is the main theorem: 

Theorem 1. Measure Pq, is continuous in a, for a G (0,1]. We precisely mean that 
Pa converges weakly to Pi as a f 1 and Pq, converges weakly to Pq^ as a ^ ao for any 
0^0 £ ( 0 , 1 ). 
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The notion of weak convergence in Theorem [T] is the classical one [1]. Theorem [T] gives 
us a new way of thinking about the stochastic heat equation with white noise. Instead 
of studying the solution to ([5]) we can study the solution to ([T]) for a 1. Also note, 
that the noise with Riesz kernel spatial covariance produces noise which is less regular. 
We like to think that this ‘roughness’ better captures properties of the stochastic heat 
equation with white noise. 

Before we begin the proof of Theorem [T], let us recall that [mild] solutions to ([T]) and 
(jSj) are interpreted as solutions of the following integral equations [6] 


Ua,t{y) 

My) 


{ua,o*pt){y) - 
{uo*Pt){y) + 



- y)a{uaM))Va{ds, dx), 
y)a{us{x))r]{ds,dx), 


where pt is the heat kernel 


Pt{x) 


1 


^/2TTKt 


exp 



2 Proof of Theorem [T] 

We will only show the first part of the theorem, converges weakly to Pi as a ^ 1- This 
is the worst case scenario. The second statement of Theorem [T] follows almost directly 
from the proof in this section. 

The proof of the upcoming Theorem [2] uses coupling, which allows us to put both 
noises rja and rj on the same probability space. This idea was introduced in [5] and lets 
us write our noise rja, for every a G (0,1) in terms of one white noise rj with covariance 

E [r]{t, x)t]{s, y)] = 6{t - s)6{x - y). 

The idea of coupling, or smoothing the noise in the spatial variable is not new. Authors 
in [3] smoothed the noise in the spatial variable by an inhnitely differentiable function 
with compact support. They have showed that this kind of smoothing converges to the 
heat equation with white noise as our smoothing function converges to 6 distribution. 
By coupling we mean that the martingale measure pa will be defined as 

PaM,t]xA)= [ [ {lA*ha,){x)p{ds,dx) 

Jo Jm. 
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and fa in ([2]) will be of the form 


fa{x) = {ha * ha){x), 

where convolution is understood as convolution of generalized functions. Function ha 
must be 


ha{x) = ci^5fi+a(a;) = gi-a (x), 

since 

9i-a{0 = 91^(0 ■91^{0- 

One might notice that ha ^ but r]a is a well dehned martingale measure, we refer 

the reader to [S] for more details. Before we state the main theorem of this section, let 
us state a technical lemma. 

Lemma 1 . | 2 [ 3 . 478 ] The following equality holds for s > 0 and 7 G [ 0 , 1 ) 

/ 1 \ -(7-i)/2 

xre-^^-^-^dx = r(-7/2 +1/2). (7) 

In the rest of this section, we will prove the following main theorem 
Theorem 2 . For every k >2 we can find 7 such that 

\imAf^^k{ua — u) = 0, 

where Af^^k is the following norm m 



Vj,fc(«) ^ sup sup (e ^‘||u(||l*(p)) , fc > 2. 

te[o,T] xeK 

Take the constant T that appears in Theorem [T] and the dehnition of the norm Af-y^k 
to be hxed throughout the whole proof. We will start our proof with Picard iterations 
for both noises rja and rj 


<r\y) 


{uo*Pt){y) + 
{uo*Pt){y) + 




y)(^{u^J'^){x)v{ds,dx) 
y)(^{u^alix))Va{ds, dx), 
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which is equivalent to the following, thanks to 0 Sec. 3.2]: 


= {uQ*pt){y)+ / {pt-s{--y)(T{u^J^\-)) *6){x)r]{ds,dx) 


Ua,t^\y) = {uQ*pt){y) + {pt-s{- - y)(y{u^a!s{-)) * K){x)p{ds, dx). 



(n)( 


First, let us estimate the L^{P) norm of the difference of Picard iterates 

nr'\y), 


E 

= E 


ui:^\y)-ur^>{y) 


(n+l) , 


( 8 ) 



{{ptsi- - * K){x) - {pt-s{- - y)(^{u^r^) * 5)(a;)) v{ds, dx) 


Adding and subtracting the following term, {pt-s{- — y)cr{u^I^'^) *ha){x) inside the integral 
and using the inequality \a — < 2^|a|^ + 2^|6|^ yields 


E 


Ad’(!/)-«r *’(!/) 


(n+l), 


< 2^E 



2''E 



{pt-s{- - y){a{u^a^s) - * h^){x)p{ds,dx) 

{pt-s{- - * {ha - 6)){x)p{ds,dx) 


The next series of steps will be used multiple times throughout this work. First we will 
use Burkholder-Davis-Gundy inequality and Minkowski integral inequality. Burkholder- 
Davis-Gundy (BDG) inequality (see for example [lOl Thm. B.lj) states that for any 
continuous martingale Mt and k >2 we have ||Mi||| < 4/c||(M)f ||fc/ 2 , where {M)t de¬ 
notes the quadratic variation of M. Applying this inequality and evaluating the quadratic 
variation term [6], Thm. 5.26] on both terms gives us 
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E 


< const • E ^ _ 

'[0,t]xR2 


,(”+1)/ 


\ fc/2' 


Pt-s{x - y)\ 3 f\x, z)fa{x - z)pt-s{z - y)dsdxdzj 


+ const • E 


'[0,t]xK2 


Pt-s{x - y){z)a{uf\x)) (/„ - 2K + 5) (x - z) 

\ k/2- 

Pt-s{z — y)o'{u^^\z))dsdxdz j 


where 


X>^:\x,z) = (a«i) - a(«f^))(x)(a«i) - a{uf^)){z). 


,W'| _ 


W'l _ 


Minkowski integral inequality states that ||/ fdp\\k < f ||/|U<^/^ for any a-finite measure 
p on M™,???. G N and jointly measurable positive function /. We use this inequality on 
the first term in to order to obtain 


E 


< const • I / 

\i[0,t]xR2 




Pt-s{x - y)vf\x, z)fa{x - z)pt-s{z - y)dsdxdz ] 


\ k/2 




+ const ■ E 


'[0,t]xR2 


Pt-s{x - y){z)a{u^J^\x)) {fa - 2h„ + 5) (x - z) 

\ k/2- 

Pt-s{z — y)(T{uf'\z))dsdxdz j 


where denotes 


vP(x, z) = E |(ff(uW) - KaK'M) - <t(uP))(z 


\k/X 


2/k 


Ultimately, we would like to show that Ua is close to m as a t" 1- The original estimate 
for Picard iterates splits into two terms 21 and 23 where 
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21 = ( / Pt-s{x - y)vf\x, z)fa{x - z)pt-s{z - y)dsdxdz ) 

'[0,t]xK2 / 


\ k/2 


® = E 


'[0,t]xM2 


Pt-s{,x - y){z)a{uf\x)) (/„ - + 5) (x - z) 


\ fc/2' 

Pt-s{z — y)a{u^^\z))dsdxdz j 


For term 21 we use Cauchy-Schwarz inequality and take supremum over the term 
involving expectation, which yields 


2t < 


sup E 




2/k 


Pt-s{x-y)fa{x-z)pt-s{z-y)dsdxdz 


The following identity holds 



(p{x)fa{x-y)ip{y)dxdy= / fa{x){ip * (:p){x)dx = / gi-a{^)\J^ip{^)\^d^, (9) 


for any p from Schwartz space iS(M) of rapidly decreasing test functions. This is a 
consequence of elementary properties of Fourier transform. We can further rewrite 21 
using identity ([9]) and the assumption that a is Lipschitz continuous with Lipschitz 
constant K, that is \cr{x) — ct(|/)| < K\x — y\ 


21 < ( / sup E 
1 0 




< 


/ sup E 

'o a:6M 




fcl 2A 


fcl2/fc f n \ 

/ gi-aiO ) 

Jr j 

\kl2 

9i-a{0\Pt-s{0f dsdn . 
; / 


k/2 


fc /2 
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Multiply by term e and obtain norm in the estimate 


[ g,.^{0 \pt-s{Of dsd^ 
') xeR L Jr 


k/2 


< 


n 






k/2 


k/2 


< - u")M / . 

) Jr I^I " / 


Later on, we will see that we can make the integral on right hand side arbitrarily small. 
The estimate for IB uses a similar technique as the estimate for 21, but some extra work is 
required because of the term (/„ — 2ha + 5) inside the integral is not a positive function. 
Thanks to 0.0 Cor. 3.4] identity (E]) extends to a much broader class of functions. We 
will use this identity to bonnd term IB. Qnantity a{u'r\-))pt-s{- — u) ^ L^(M) n L^(M) 
almost surely, because 


E 




< 2K\1 + Af,,,{u^-^r)\\pUWLHR) , 


for some 7 > 1 and N'^^ 2 {udJ) is uniformly bounded for every n G N uni. A similar 
reasoning applies for H-ULqK). We can write 


E 


= E 


' Pt-s(x - y)a(ud"\x)) {fa - 2ha + 6) {x - z)pt-s{z - y)a{u^f^\z))dsdxdz ) 

[ 0 ,t]xR 2 / 

\ k/2 

{91-a - 2gi^ + 1)(0 \J^ {pt-s{- - y)(^{u^r\-))) 


\ k/2 


' [0,4] xR 




Split this integral into two parts, and use inequality |a — 
to get 



< const • ('^ + 2 )) ) 


\ k/2 

2 g(i-a )/2 + 1)(0 1 -^ {pt-s{- - y)(T{u^r\-))) {0\^ d^dsj 













where 


e = E 

D = E 


'[-1.1] 


{gi-a - 2 ^ 1 ^ + 1)(0 1 *^ {pt-s{- - y)(j{u^J^\-))) ( 0 | d^ds 


\ k/2 




[- 1 , 1 ] 


{gi-a - 2gi^ + 1)(0 1-7^ {pt-s{- - 1 /)CT(«(”)(•))) ( 0 | d^ds 


\ fc/2' 




Properties of Fourier transform and Lipschitz continuity of <t{x) give us 

\T - !/)<r(«<”l(-))) iOf < IIP.-.(- - !/)<^(d”'(-))lli.(R) 

< A'=||p,_,(- - p)(l + |«<”'(')l)lli.(.| < + 2||p,-.(- - P)«[”’(-)lli.(E,). (10) 

for the term inside of Splitting the term fIlOp inside of the integral into two yields 


(J:< E , , 

'[0,1] 2[-1,1] 
< Ca + const • E 


(^i_„ - 2gi^ + 1){^)K^{2 + 2\\pt-si- - yWr\-)\\hiR))d^ds 


k/2 


'[ 0 ,i[ 2 [- 1 , 1 ] 


{gi-a - 2^1^ + i)(0lbt-s(- - y)u‘s'\-)\\lHR)d^ds 


k/2 


where Ca denotes 

Ca = const 


'[o,t] 2 [- 1 , 1 ] 


{gi_a - + l){0d^ds 


k/2 


Term Ca can be made as small as we like, due to the dominated convergence theorem. 
We use Minkowski integral inequality and get 


<Ca + const 
<Ca + const 


/ sup E 

'[0,i] xGR 




2/k 


/ e-'=^*supE 

'[0,1] 3;GR 


/ (^l-a - 2^1^ + 


'[- 1 , 1 ] 

2/k 


{gx-a - 2gi^ + l){i)dids 


k/2 


[- 1 , 1 ] 


k/2 


<Ca + const • [ {gi_a - 2,gi^ + l){^)d^ds\ 


k/2 


From the general theory of stochastic partial differential equations im we know that the 
term is bounded uniformly for every k, for some 7 sufficiently large. The integral 
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term bounding CC can be made arbitrarily small, again from the dominated convergence 
theorem. 

All we have left to do is hnd the estimate for 2). Add and subtract the term cr(ns(s)) 
inside the Fourier transform, split into two integrals and obtain 


2) < const • sup E 

n>no 


+ const ■ E 



{pt-s{x - y) {a{uf\x)) - a{us{x)))) dxdtj 



{gi-a - 2 ^ 1 ^ + 1)(0 \J^ (Pt-s(- - y)(T(us(-)) ( 01 ^ d^ds) 


kji 


0 ^'RVhl,!] 




for large n> uq. The hrst integral can be made arbitrarily small, this is from convergence 
of Pickard’s iterations and theory of SPDEs [10]. We have used Plancherel’s theorem 
and the fact that {gi-a — “^gi^ + 1) is bounded by constants on M \ [—1,1], uniformly 
for all a G (0,1). The second term can be made arbitrarily small due to the dominated 
convergence theorem. If n is small, that is smaller than some no, then we need to pick 
appropriate a close to one so that the term 2) is arbitrarily close to zero for all n < uq. 
This again follows from the dominated convergence theorem. The term (qi_^ —2qi-a +1) 
inside 2 converges pointwise to zero. 

The estimate of is bounded by the norm of the previous Picard iterate. 

Both terms € and 2 can be made smaller than arbitrary e, thus term 53 can be made 
smaller then arbitrary e. Add up all those estimates and multiply ([HD by to get 


-k'yt 




(n+l), 




k/2 


< const ■ ( / e-^^d-s) / J + ^ ^ 

D Jm. I^I / 


where 7 > 1. We use Lemma [T] to evaluate the term inside the integral. A straightforward 
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calculation yields 


/ \ 

< const ■ U “ s)-“/"dsJ + 


< const ■ + 


k/2 


< const ■ N'^^kiu^a^ — 


l-a/2\ 


/ 1 \ fc/4 

+ e < const ■ f - J + 


e . 


We can take supremum over y eM. and t E [0, T] to get 

/ 1 \ fc/4 

< const • ( -j + e . 


/, \ fc /2 

This defines a convergent geometric series assuming that the coefficient C I-j < 1, 
where C = const. This also implies 

— ~ Q 5 

1 “ :jy72 

for some a close to 1. This concludes the proof of Theorem O 


2.1 Convergence of finite dimensional distributions 

Theorem [2] also states that the solution converges to u in L‘^{P) norm for every 
t E [0,T] and x G M. This implies weak convergence of finite dimensional distributions 
of Ua to finite dimensional distributions of u. The easiest way to see that is to show 
convergence in probability for a finite number of pairs (tj, Xj) G [0, T] x M, which implies 
weak convergence of finite dimensional distribution [H pg. 27]. The convergence in 
probability follows from Chebyshev’s inequality and convergence in T^(P) norm. 

2.2 Estimates for Kolmogorov’s continuity theorem and Tight¬ 
ness 

We will prove tightness (and thus weak convergence) from Kolmogorov’s continuity the¬ 
orem. Before we begin the proof, we will need the following two lemmas. 
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Lemma 2. Lemma 6 . 4 ] For all s > 0 and x G M 



~ Pt{y) \ dy < const • 



where the implied constant does not depend on {s,x). 
Lemma 3. For all t,e > 0 and x G M we have 


\pt+e{y) - Pt{y)\dy < const ■ ((log(t + e) - log(t)) A 1) 


Proof. Direct compntation gives ns 

\pt+e{y) -Pt{y)\dy = 

/*t+e /* / 1 „,2 



Ps{y)ds 


< 


dy = 





h + iL- I pAy)ils 




2t 2 Fk 

dt = (log(t + e) - log(t)) 


dy 


In addition, we have that 


\pt+e{y) -Pt{y)\dy < 2. 


□ 


2.2.1 Difference in the spatial variable 

Let us estimate the difference in the spatial variable 


E 


u. 


a,t 


X] 


U 




iy)\' 


= E 



{pt-s{x - z) - Pt-s{y - z)) a{ua,s{z))r]{ds, dz) 


and denote 


B{z) = {pt-s{,x - z) - pt-siy - z)) , 
A{x,y) = (T{ua,s{x))a{ua,s{y))- 
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We will proceed just as in Section |21 We use BDG, Minkowski integral inequality, 
Cauchy-Schwarz inequality and take the absolute value inside the integral and get 


E 





fa{z — w)B{z)B{w)A{x, y)dsdzdw 


fc/2' 



k/2 


k/2 


< const / sup \\a{Ua,s{x))\\k{fa *Pt-s){0) / \pt-s{x - z) - Pt-s{y - z)\dsdz 

Jo Jr 

nt p k/2 

< const(l/ {fa*Pt-s){0) / \pt-s{x - z) - pt^siv - z) \ dzds 

Jo Jr 

A f\\x-v\\ ^ 

< COnst(l -h J {fa*Pt-s)iO) i - - A 1 


y/kj: 


where the last inequality is due to Lemma [2j We also used the fact that 

/ fa{z - w) \pt-s{x - w) -Pt-s{y -w)\dw <2 {fa * Pt-s) (0). 


The inequality r A 1 < for a G (0,1/2) gives us 


E 



Ua,t{y)\ 


k 


< const(l + Jfj,k{Ua)’")\\x - y\\ 


ak 


{fa*Pt-s){0) ■ {t- s) ‘"ds 


k/2 


( 11 ) 


^0 

It remains to show that the integral on the right hand side is bounded for all a G 
(ctoj 1), «o > 0. To show this, we will need an explicit form of fa- The result is stated in 
the next lemma. 

Lemma 4. For every 1 > a > Oq > 0 we have 

fa*Ps{0) < const • 

where the constant depends only on our choice of Oq. 
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Proof. By direct computation and ([7]) we get 

■ 1 


(/a *Ps)(0) = Ci_, 


Ct / I IQ: 

P \x\ 


Ps{x)dx = 2 


sin I ) r (a) fl-a\ 

- 2-PT i s-“/27r-V2 


(27r; 


a 


V 2 


The boundedness of constant 


sin 


(1—a)7r 


r(a) 


2-t./2r 


1 — a 


(27r)“ ‘ V 2 

can be concluded from Euler’s reflection formula. 


71 


- 1/2 


□ 


Because of the lemma above, the integral on the right hand side of (ITT]) is hnite as 
long as a G (0,1/2). 

2.2.2 Difference in the time variable 

The difference in the time variable is going to be, for 5 > 0 


E 


\Uc,t+5{x) - Uc,t{x)\ 
= const ■ E 


t p k' 

{pt+&-s{x - z) - Pt-s{x - z)) a{ua,s{z))pids, dz) 



+ const • E 


pt-\-5 p ^ 

/ Pt+s-s{x - z)a{ua,s{z))p{ds,dz) 


Let us estimate the second integral, we can use the same technique as in the case of the 
spatial variable and write 


E 


r*t+5 p k 

/ Pt+5-six - z)a{ua,siz))p{ds,dz) 


rt+S 

< const I / sup E 


<7 


'Ua,si^))\ 

•t+5 


lei 


\pt+s-s{0f dfds 


k/2 


/ rt+a 

< const(l + i {t + 6 — j 

< COnst(l + N'j,k{Ua)’^) < COnst(l + Af'y,k{Ua)^) 
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The estimate for the second integral will be 


E 



{pt+s-s{x - z)- Pt-s{x - z)) a{uo,^s{z))p{ds, dz) 

2/A: r 


ft 

r 1 

/ sup E 

\a{Ua,six))\’' 

Jo xGR 



\k/2 

V 


< const(l + ( / {fa*Pt-s){0) / \Ps+6{z) - Ps{z)\ dzds 

\Jo Jr 

/ \ 

< const(l + A/'^,a;(m«)^) {j (log(s + (5) - log(s)) dsj 


k/2 


< const(l + ^4\/5atan + 2\/tlog(l + (5/t) 


k/2 


( 12 ) 


by using a similar technique as in the case for the spatial variable and Lemma [3l The 
inequality log(l + C) < vT for all C > 0 gives us 


E 


t p ^ 

{pt+5-s{x -z)- Pt-s{x - z)) a{ua,s{z))'n{ds, dz) 



< const + [A/"^,A:(Wo)]^j (13) 

We can combine both estimates flT^ and flT^ to finally get 


E 


u, 


,,t+six) - Ua,tix)\^ = const ■ (^1 + [A4,fc(Ma)]'') • 


2.3 Kolmogorov’s continuity theorem and tightness 

Let us mention that Nk,^{ua) is finite O Prop. 9.1] for every choice of a G (0,1). We 
also know that Ua is continuous in norm thanks to Theorem [3l Those two facts 
together imply that M^^kiuo) is uniformly bounded for a G (ao, 1), ao > 0. We would like 
to emphasize that we are not creating a circular argument, since the proof of Theorem [3] 
solely depends on the proof of Theorem [2l We have that for every 1 > a > ao > 0 and 
(s, x), (f, y) from [0, T] x [—N, iV] C Mq x M the following holds for fc > 2 


E 


\^a,s(^x) Ua^t(^y)\ 


< const lx 


+ const \t — , 
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where a e (0,1/2), thanks to our estimates from sections r 2 . 2 . 1 l and l 2 . 2 . 2 [ Denote pit, x) = 
|x|“ + , then Kolmogorov’s continuity theorem states that there is a modihcation of 

Ua,s{x) such that (see for example Uni pg. 113]) 


E 


sup 

{s,x),{t,y)eD 


Ua,s{x) - Uq^tjy) 
p{s -t,x- y)i 


k 


< A < +00 


(14) 


for every a G (ao, 1) and q G (0 ,1 — H/k) where H = l/a + 4. By Chebyshev’s inequality 
and flT^ . we can write 


which implies 


P I sup |Ma,s(a;) - Ua,t{y)\ > ^ 

{s,x),{t,y)eD 

\p{s-t,x-y)<5 


< 4 ' 5 *’ . 


lim sup P 

(S^-O Qg(c(,^i) 


{ sup \Ua,s{x) - Ua,t{y)\ > e 
{s,x),{t,y)£D 
p{s-t,x-y)<S 


0 , 


for every e > 0. We can conclude [I3l Thm. 2] that the solution converges weakly to 
u in C. 


3 Continuity in a for a G (0,1) 

Our proof in section [2] also implies continuity in a for a G (0,1). We will only comment 
on how the proof would change in section [ 2 l 

Results in sections 12.2.2112.2.11 and 12.31 hold without change. The piece of proof that 
needs to be slightly modified is the proof of Theorem |2l Recreating Theorem |2] for the 
new setting would not give us any new technique or insight. Let us simply state the new 
Theorem. 

Theorem 3. For every k >2 and oq ^ (0,1) we can find 7 such that 

lim Af^^kiua - M«o) = 0 - 

The proof of Theorem 3 follows the same general direction of the proof of Theo¬ 
rem [2] with the following changes. We need to replace ufix) with Uao,t{.x) and change 

{fa — 2ha + 5) in the estimate for to — 2f a+ao + fa^^ and change (^gi-a — 2,gi^ + 1 
lo {gi-a — ‘^gi_st+^ + gi-ao^ ■ 
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